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The three neutrino model has 9 physical parameters, 3 neutrino masses, 3 mixing angles and 3 CP violating 
phases. Among them, neutrino oscillation experiments can probe 6 parameters: 2 mass squared differences, 3 
mixing angles, and 1 CP phase. The experiments performed so far determined the magnitudes of the two mass 
squared differences, the sign of the smaller mass squared difference, the magnitudes of two of the three mixing 
angles, and the upper bound on the third mixing angle. The sign of the larger mass squared difference (the 
neutrino mass hierarchy pattern), the magnitude of the third mixing angle and the CP violating phase, and a 
two-fold ambiguity in the mixing angle that dictates the atmospheric neutrino oscillation should be determined 
by future oscillation experiments. In this talk, I introduce a few ideas of future long baseline neutrino oscillation 
experiments which make use of the super neutrino beams from J-PARC (Japan Proton Accelerator Research 
Complex) in Tokai village. We examine the potential of HyperKamiokande (HK), the proposed 1 Mega-ton 
water Cerenkov detector, and then study the fate and possible detection of the off-axis beam from J-PARC in 
Korea, which is available free throughout the period of the T2K (Tokai-to-SuperKamiokande) and the possible 
T-to-HK projects. Although the CP violating phase can be measured accurately by studying — > v e and 
F M — > 77 e oscillations at HK, there appear multiple solution ambiguities which can be solved only by determining 
the neutrino mass hierarchy and the twofold ambiguity in the mixing angle. We show that very long baseline 
experiments with higher energy beams from J-PARC and a possible huge Water Cerenkov Calorimeter detector 
proposed in Beijing can resolve the neutrino mass hierarchy. If such a detector can be built in China, future 
experiments with a muon storage ring neutrino factory at J-PARC will be able to lift all the degeneracies in the 
three neutrino model parameters. 



1. Introduction 

All the experimental observations on neutrino 
physics, with a notable exception of the LSND 
experiment [1], can be accommodated well in the 
three neutrino model. The atmospheric neutrino 
oscillation first observed by Super-Kamiokande in 
1998 [2] and the observed deficit of the muon neu- 
trino flux in the K2K experiment [3] can be in- 
terpreted as — > v T oscillation with the 90% CL 
allowed ranges [4] 

1.9 x 10~ 3 < Sm 2 (eV 2 ) < 3.0 x 10~ 3 , 

ATM v ' 7 

sin 2 20 ATM >0.90. (1) 

The observed deficits of the solar neutrinos [5] 
and the reactor anti-neutrinos by the KamLand 
experiment can be interpreted as v e — ► or v T 
oscillation with [6] 

Ko L = 8-218:5 x 10- 5 eV 2 , 



t an20 so L -°- 4O -ao 9 7- ( 2 ) 

No observation of the deficit of the reactor V e flux 
constrains the third mixing angle [7] 

sin 2 20 RCT < 0.20 for 5m 2 = 2.0 x 10" 3 eV 2 , 
sin 2 26» RCT < 0.16 for 6m 2 = 2.5 x 10" 3 eV 2 ,(3) 
sin 2 26» RCT < 0.14 for Sm 2 = 3.0 x 10" 3 eV 2 , 

at the 90% CL. Since the magnitude of the mass 
squared difference of order 1 eV 2 associated with 
the — ► v e oscillation observed by the LSND ex- 
periment is much bigger than those of the above 
differences, once their finding is confirmed by the 
MiniBoone experiment [8], we need at least one 
more neutrinos. Nevertheless, I would like to dis- 
cuss physics prospects of future neutrino oscilla- 
tion experiments in the framework of the three 
neutrino model, mainly because of my prejudice 
that I find no compelling theoretical reason for 
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the existence of more than three light neutrinos. 
Therefore, once the LSND observation is con- 
firmed by MiniBoone, we should revise all our 
thoughts from scratch. 

In the three neutrino model, the three weak in- 
teraction eigenstates v a (a = e, [i, r) that appear 
in the universal charged current interactions and 
the three mass eigenstates Vi (i — 1,2,3) are re- 
lated by the 3 x 3 MNS (Maki-Nakagawa-Sakata) 
matrix [9] 
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The MNS matrix can be expressed as [5] 

V uss =UV = U&mg(e^l\e^l\ 1) (5) 

with the two Majorana phases a.\ and a 2 . The 
3x3 unitary matrix U can be parameterized in 
terms of the three angles and one CP violating 
phase just like the CKM matrix: 



u = 23 P s 13 P}0 12 

( u el u e2 u e3 

U Ufa U ^3 
y U T i U T 2 U t3 



(6) 



where Oij are the orthogonal rotation matrix in 
the ij plane, and Ps — diag(l, 1, e lS ) gives the CP 
violating phase of the MNS matrix 5 — <5 MNS . Be- 
cause the present neutrino oscillation experiments 
constrain directly the three elements, 



u e3 


-iS 

= si 3 e 


(7a) 


u e2 


= S12C13 , 


(7b) 


U^ 3 


= S23C13 , 


(7c) 



we find it most convenient to adopt the conven- 
tion [10] where these three matrix elements in the 
upper-right corner of the U matrix are chosen as 
the independent parameters. The phase conven- 
tion of Refs. [10,11] 



U e3 = \U e3 \e~ lS 

Uel > U e2 > , 

u„ 3 ,u t3 > 0, 



(8a) 
(8b) 
(8c) 
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Figure 1. Schematic view of the four-types of 
neutrino-mass hierarchy. 



is equivalent to the choice < #12 < 7r/4 and 
< #13, #23 < 7r/2. In our convention v\ is de- 
fined by the condition U e \ > U e2 , and hence both 
mi < m 2 and mi > m 2 are allowed, and they 
should be determined by experiments. *Four pos- 
sible patterns of the neutrino mass hierarchy are 
then defined as [11] 



(Sm 2 12 ,Sm 2 13 ) 
= ( Sm L^ 



(9) 



5m 
5m 



) (I: normal hierarchy) 



2 

ATM 

Lm) (II) 

2 

ATM 



= ( 5m 2 OL , -(Sm^ TM ) (III: inverted hierarchy) 

= (-KoL."K™) (IV) 

where Smfj = m 2 — mf . See Fig. 1. The solar 
neutrino experiments determine 
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= m 2 -m 1 = <5m soL > 0, 



(10) 



and hence only the hierarchy patterns I (normal) 
or III (inverted) are yet to be determined by fu- 
ture experiments. 

We keep the four hierarchy cases in the follow- 
ing discussions because we find the following the- 
orem [11,12] useful: 

P-t: = P?, V ^.,,,. (normal hierarchy) , (11a) 
(inverted hierarchy) (lib) 



pll 



"This is equivalent to an alternative definition of u\ where 
mi < ni2 is assumed and the ordering U e i > U e 2 or V e \ < 
U e 2 is to be determined by experiments. 
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The theorem is valid in the presence of an arbi- 
trary matter profile [12]. Because of the theo- 
rem (11), we do not give expressions for the anti- 
neutrino oscillation probabilities explicitly. For 
instance, Pp ^,v c in the normal hierarchy is ob- 
tained from the expression for P v ^, Ve simply be 
reversing the sign of the factors A i2 and A13. 

The probability of finding the flavor-eigenstate 
(i from the original flavor-eigenstate a at the base- 
line length L in the vacuum is given by 



Vet— 

3 



(12) 
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6m 
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5m% (eV 2 ) 
lJ L~ 2.534 „/ £(km). (13) 



E v (GeV) 



It is clear from this expression that the neutrino 
oscillation experiments cannot measure the two 
Majorana phases ol\ and 0:2 in cq. (5), and that 
only the two mass squared differences out of the 
three neutrino masses can be probed by those ex- 
periments. 

It is useful to give expressions of the oscilla- 
tion probabilities when either A13 or A12 is order 
unity. For instance when | A.13 1 <~ 0(1), we can 
expand the expressions in terms of | A.12 1 <C 1 and 
find 



p 

+ 21^2 
Pu.^u. = 1 



l-4|^ 3 | 2 (l-|t^3 
M 3| sin A13 • A12 



! )sin 2 



2 \U, 12 



0((A 12 ) 2 ),(14a) 

A13 

2 

+2|[/ e2 | 2 |[/ e3 | 2 sin A13 • A 12 + 0((A 12 ) 2 ) , (14b) 
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l^esHsin 2 



//3 



2\tt |2 • 2 ^13 



+2Re{U^U: 2 U e3 li; 3 } sin A 13 • A 12 



+4J s sin 2 ^.A 12 + 0((A 12 ) 2 ). 



(14c) 



The expression eq. (14a) is relevant for the atmo- 
spheric neutrino oscillation and the K2K experi- 



ment, and wc may identify 



sin 2 29 
dm 



ATM 
2 

ATM 



4|C/ M3 | 2 (1-|^3| 2 ), (15a) 
\Smj 3 \ = |m§ - m\\ . (15b) 



The expression eq. (14b) is relevant for the reac- 
tor anti-neutrino oscillation experiments after the 
replacements in eq. (11) is made: 



sin 2 26 
5m 



RCT 
2 



4|t/ e3 | 2 (l-|E/ e3 | 2 ), 



1 5m 2 



131 



5m 2 



(16a) 
(16b) 



Only the magnitude of 5m 2 3 is constrained be- 
cause, so far no effects of order A i2 have been 
positively identified. Finally, cq. (14c) is relevant 
for the — > v e appearance experiments such as 
T2K. It is important to note that the — > v e 
oscillation probability is proportional to 



4 sin' 



#atm sin 



411^3 11^ 



(17) 



and that the CP violating effect appears in the 
A12 order. The Jarlskog parameter of the MNS 
matrix is defined as 



J 



= lm{U ai U* pi U 0j U* aj ) 



i-\u, 



Im(J7 B3 ) 



(■:.', 



= C 2 3 Ci2C 2 3Sl2S23Sl3Sin(5 Iv 



(18a) 
(18b) 
(18c) 



where the orderings (a, (3) = (e, fj,), (fi, r), (r, e) 
and (i, j) = (1, 2), (2, 3), (3, 1) should be taken in 
the definition. The expression (18b) follows from 
the convention (8) [10], and eq. (18c) follows from 
the parameterization (6) [5]. 

Comparison of the expansions in eq. (12) and 
the theorem (11) tells that the anti- neutrino os- 
cillation probabilities have the same expressions 
as those for neutrinos up to terms of odcr A 12 , 
except for the term proportional to J MNS - We 
find 



P — P- - 

= 8J MNS sin 2 ^H.A 12 + 0((A 12 ) 2 ), (19) 

which is valid for both normal and inverted hi- 
erarchies. Although the above expression is valid 
only when we can neglect the earth matter effects 
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of the oscillation probabilities, it gives a good ze- 
roth order approximation to the proposed T2K 
and T-to-HK experiments, which plan to use sub- 
GeV neutrino beams. 

When the magnitude of A i2 is of the order of 
unity, |A 13 | - Sm 2 ATM /5m 2 SOL ~ 30 > 1 applies, 
and we have an alternative expansion of eq. (12). 
The most relevant one is 

P,.-,. =l-4|l/ cl p|<7 e3 |W^i 

where SE^/E^ is the experimental resolution of 
the neutrino energy. The expression (20) is 
relevant for the solar neutrino oscillation and 
the KamLand reactor neutrino oscillation experi- 
ments. By neglecting the small term proportional 
to sin 2 2# RCT , we can identify 



sin 2 26. 
5m 



SOL 
2 

SOL 



4|C/ e2 | 2 |C/ el | 



\5m 



12 



m 9 



m. 



(21a) 
(21b) 



The sign of the mass squared difference cannot 
be determined from the above expression, which 
applies for the KamLand experiments. 

From the above arguments, we can parameter- 
ize the complete MNS matrix elements by using 
the observed constraints on the three mixing an- 
gles: 



1 \ 1 



u, 



e2 



2 _ 1 - \U e s\ 2 V / (l-|t/ e3 | 2 ) 2 -sin 2 20 soL 



1 ± ^jl - sin 2 29 ~ 



(22) 



Since the magnitudes of the three matrix ele- 
ments at the corner of the MNS matrix are de- 
termined as above, the full MNS matrix can be 
determined for a given CP phase (5 MNS [10]. Al- 
though the above identifications are valid only ap- 
proximately, they suffice for our discussions about 
physics prospects of future neutrino oscillation 
experiments. 



2. Matter effects in the neutrino oscillation 

The neutrino mass hierarchy pattern between 
I and II (III and IV) has been distinguished by 
the matter effects of the observed solar neutrino 
oscillation. 

The oscillation probabilities in the matter has 
the same form as in eqs. (12), (13) and (14), where 
the mass-squared differences and the MNS matrix 
elements are replaced by those in the matter: 



A,- 



•A 



(23) 



if the variation of the matter density remains 
small along the base-line. 

Neutrino-flavor oscillation inside of the matter 
is governed by the equation [13,14] 
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(24a) 



where the Hamiltonian (in units of 2E V ) in the 
matter is 



H = U diag(0, 5m 2 12 , 5m 2 13 ) C/ f 
+diag(a,0,0) 

/ Ai \ 



= U 



A 2 






A 3 



(24b) 
(24c) 



Here a is the matter effect term, 
a = 2V2G F n e E u 

= 7.56 x 10- 5 (cV 2 ) ( — f - 
\ g/cm 



ik) << 25 > 



where n e is the electron density, Gf is the Fermi 
constant, and p is the matter density. 

The Hamiltonian H in the matter is diago- 
nalized by the MNS matrix in the matter U. 
The neutrino-flavor oscillation probabilities in the 
matter take the same form as those in the vacuum 
by making the replacement (23) where 
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(26) 
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Figure 2. The magnitudes of the effective mass 
squared differences, |#mf 3 | (l^m^l in the hierar- 
chy II, IV cases) and |<5m 2 2 |, as functions of the 
neutrino energy times the effective matter density 
Pes (g/cm 3 )/3. In the earth crust, p e g ~ 3, and 
at the center of the sun, p e s ~ 200. 



In Fig. 2, we show the ^-dependence of the mag- 
nitudes of the effective mass-squared differences, 
|(5m^ 3 | and |<5m 2 2 | for the 4 mass hierarchy cases 
of Fig. I. In case of the hierarchy II and IV, 
we show |<5to 23 | rather than |£m 2 3 |, because the 
leading oscillation term is dictated by the former 
difference. The curves are obtained for 



|<5m 2 3 | 
\5m 2 12 \ 
tan 2 0„ 



sin 



2.5 x 10~ 3 eV 2 . 



8.3 x 10" 
= 0.4, 
= 0. 



(27a) 
(27b) 
(27c) 
(27d) 



It is worth noting that the matter effects depend 
only on the above four parameters, i.e., the two 
mass squared differences and the two angles, 0i 2 



and 0i3, because of the identity [15] 
H = 2S P 5 H'Pl0l 3 , 

where 



(28a) 



H' = 13 12 diag(0, 5m 2 12 , 5m 2 13 ) O t 12 O t 13 

+ diag(a,0,0). (28b) 

It is clear that the above reduced Hamiltonian H' 
can be diagonalized as 



H' = Odiag(Ai,A 2 ,A 3 )0 3 



(29a) 



with the same eigen values of H in eq. (24b). The 
MNS matrix in the matter is then expressed as* 



U = 23 P s O. 



(29b) 



In the limit of sin 2 rct = (s 13 ) 2 = in cq. (27), 
the matrix 0\ 3 = 1, and the eigen value prob- 
lem is reduced to the 2x2 matrix problem. The 
eigenvalues are 



A± 
A 3 



8m\ 2 + a ± {5m\ 2 + a) 2 — Aac\ 2 5m\ 



12 



5m\ 3 , 



(30) 



and the results are shown in Fig. 2 for all the four 
hierarchy cases of Fig. 1. Note that 

(Ai,A 2 ) = (A_,A + ) hierarchy I, III, (31a) 
= (A+, A_) hierarchy II, VI , (31b) 

so that Ai — > and A2 — Ai = 5m\ 2 — > 5m\ 2 in 
the a — > limit. Although all the curves in Fig. 2 
are obtained in the limit of sin 2 RCT = 0, we find 
that the qualitative features of the curves remain 
the same even when sin 2 20 RCT w 0.1, except near 
the 10 GeV region. 

The horizontal scale of Fig. 2 is the neutrino 
energy times the effective matter density in units 
of 3 g/ (cm) 3 , and the vertical scale gives the mag- 
nitudes of the effective mass squared differences 
in the matter. The typical energy of the reac- 
tor anti-neutrinos is about a few MeV, or a few 



^Note that once the expression (29b) is used to calculate 
the MNS matrix in the matter, the phase convention of 
eq. (8) is no more valid. 
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x 10~ 3 GeV, and hence there is no significant dif- 
ference between the hierarchy I and II (or between 
III and IV) in the CHOOZ or KamLand exper- 
iments. On the other hand, the typical Boron 
( 8 B) neutrino flux energy of 7 MeV with the ef- 
fective matter density of p cS < 200 g/(cm) 3 in 
the sun extends up to 450 MeV in the horizontal 
scale, where the matter effects clearly distinguish 
between the neutrino mass hierarchies, and the 
energy dependence of the observed solar neutrino 
flux clearly favors the hierarchy I or III, over II or 
IV. In fact, in case of the hierarchy I or III, the 
effective mass squared difference 

\5m\ 2 \ = |A 2 - Ai| (32) 
= \J {a - cos 2#i2 8m\ 2 ) 2 + (cos 29i 2 8m\ 2 ) 2 

has the minimum at 

a = cos 2#i2 5m\ 2 (33) 

where the effective mixing angle 9\ 2 crosses 7r/4. 
On the other hand, in case of the hierarchy II or 
IV, since 5m\ 2 < in eq. (32), there is no mini- 
mum, and the effective mixing angle decreases at 
high energy or density. The difference is signifi- 
cant enough to affect the observed energy spec- 
trum of the solar neutrino flux on the earth, and 
the neutrino mass hierarchy I or III has been cho- 
sen against II or IV. 

Between the hierarchy I (normal) and III (in- 
verted), the difference in the matter effects re- 
mains small for the effective mass squared dif- 
ferences in Fig. 2. The ongoing K2K experi- 
ment with the baseline length of 250 km and the 
planned long baseline neutrino oscillation experi- 
ments T2K, from J-PARC at Tokai to Kamioka, 
whose baseline length is about 300 km, adopt 
rather low energy neutrino beam below 1 GeV, 
where the matter effect is rather small. We will 
find below that the remaining ambiguity between 
the neutrino mass hierarchy I (normal) and III 
(inverted) can be distinguished by using the earth 
matter effects at higher energy very long baseline 
neutrino oscillation experiments. 



3. T2K and T-to-HK projects 

The proposed long-baseline (LBL) neutrino- 
oscillation experiments [16,17] will achieve the 
precision measurements of <5ra 2 TM and sin 2 2# ATM 
by using conventional neutrino beams, which are 
made from decays of tt and K that are produced 
by high-energy proton beams. In particular, the 
T2K experiment will be able to detect — > v e 
appearance, if the magnitude of the oscillation 
eq. (17) is not too small. 

The T2K neutrino beam line is now under 
construction at J-PARC (Japan Proton Accel- 
erator Research Complex) in Tokai village [18], 
whose 50 GeV high intensity proton accelera- 
tor with 0.75 MW will deliver 10 21 POT (pro- 
ton on target) per year. This is about a factor 
of 50 more intensive than the KEK proton syn- 
chrotron for the present K2K experiment. Since 
Tokai village is about 50 km east from KEK, 
the baseline length of the T2K experiment will 
be L — 295 km, as compared to L = 250 km 
for K2K: see Fig. 3. They plan to make use of 
the K2K proven technique to reconstruct neu- 
trino energy at Super-Kamiokande (SK) by se- 
lecting quasi-elastic charged current events. The 
expected accuracy is 

A(|<5m 2 3 |) = 10~ 4 eV 2 , (34a) 
A(sin 2 20 ATM ) = 0.01, (34b) 

from the precision measurement of the deficit of 
the flux. Because SK is capable of detecting 
v e charged current events clearly, they will have a 
good sensitivity to the — ■* v e appearance mea- 
surement if 

sin 2 29 r , > 0.006. (35) 

RCT ~ I 

Even with the intensive proton beam from J- 
PARC, it is not an easy task to measure the oscil- 
lation probabilities for anti-neutrinos because the 
charged current cross section of anti-neutrinos is 
a factor of 3 smaller than the neutrino cross sec- 
tions. A factor of 5 upgrade in the power of the 
proton beam to 4MW, and a factor of 50 larger 
detector, 1 Mton water Cerenkov detector Hyper- 
Kamiokande (HK) are therefore envisaged [19]. 

In Rcf.[20], we studied the possibility of mea- 
suring 5 in the future T-to-HK experiment. 



Figure 3. The long baseline neutrino oscilla- 
tion projects, K2K (KEK to Super-Kamiokandc; 
250km) and T2K (J-PARC at Tokai to Super- 
Kamiokande; 295km). 



There the measurement of the difference (19) by 
making use of both the neutrino [y^ enriched) 
beam and the anti-neutrino (v^ enriched) beam 
allows us to study CP violation in the lepton sec- 
tor. We find that the advantage of the T-to-HK 
experiment over various proposed neutrino oscil- 
lation experiments lies in the fact that the water 
Cerenkov detector is capable of measuring both 
H and e charged current reactions, and that the 
use of low energy beam (of order 1 GeV or less) 
at relatively short distances (295 km) makes the 
matter effects small. Because of the above advan- 
tages, the experiment can effectively measure the 
difference (19), and hence have high sensitivity 
to the CP violating leptonic Jarlskog parameter 
J t of eq. (18), or sin 5 » For instance, it is 

MNS ^ \ n MNS ' 

relatively easy to distinguish between J MNS = 90° 
and S c = 270°. On the other hand, we found 

MNS " 

that it is a non-trivial task for the experiment to 
distinguish between 



most notably, between n = 0° and <5 t „ = 

J ' MNS MNS 

180°. In order to distinguish between the two 
cases in eq. (36), which give the same sin <5 MNS , we 



Figure 4. The neutrino beam flux from J-PARC 
at L=295km for different off-axis angles, 3, 2.5, 
2, 1.5, 1, 0.5 and degree (on axis). 



should measure cos S, ,„ TD as well. In the v u — > v e 

MNS P c 

transition probability expression of eq. (14c), the 
relevant term is 

2Re[U^ 2 U: 2 U e3 U; 3 ] sin A 13 • A 12 , (37) 

whose magnitude is small at around Ai 3 <~ n 
where the leading sin 2 (A13/2) term takes its max- 
imum value of unity. In order to distinguish this 
sub-leading term from the leading one, we need to 
make experiments at different neutrino energies. 
By choosing the neutrino narrow band beam at 
two different energies, we could show in ref. [20] 
that the degeneracy between 5 MNS = 0° and 180° 
can be resolved. 

After our work [20] was completed in 2002, the 
T2K collaboration decided to adopt the off-axis 
beam (OAB) rather than the narrow band beam 
(NBB) which was used in our analysis. The OAB 
flux has a sharp peak at low energies and has 
relatively long high-energy tail. The flux inten- 
sity of the OAB around the sharp peak is typ- 
ically stronger than that of the NBB, while the 
hard high energy tail may be considered as a dis- 
advantage. We show the OAB flux profile from 
the J-PARC 50 GeV proton accelerator in Fig. 4. 
Fluxes at off-axis angles between 0° (on-axis) and 
3° are shown for later use. For the T2K and T- 
to-HK experiments, the neutrino beam line at J- 
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PARC has been set up such that the OAB beam 
between 2° and 3° can be sent to SK and HK. 
The beam orientation has been so chosen that 
once HK is constructed at its planned location, 
exactly the same OAB will hit both detectors. 
We may read off from the figure that the peak 
flux energy is about 0.8, 0.7, and 0.6 GeV for 2°, 
2.5°, and 3°, respectively. The larger the OAB 
angle, the lower the peak energy is, and sharper 
the cut-off for the high-energy tail. 

By repeating the analysis of ref. [20] for the 
experimental setup of 

0.8 Mton- year for ^OAB(2°), (38a) 
2.0 Mton- year for ^OAB(3°), (38b) 
3.2 Mton • year for 77 Al OAB(2°) , (38c) 

we find the results of Fig. 5. The flux ratios 
are so chosen that roughly the same number of 
charged current events are expected, about 10 4 
/i-like events and about a few hundreds to a thou- 
sand e-like events, in the range of the parameter 
space we studied. If a factor of 5 more inten- 
sive beam is available at J-PARC by the time we 
will have HK, the same quality results may be 
obtained in less than one and a half year time. 

The contours of Fig. 5 show the 1-, 2-, and 3-cr 
allowed region in the plane of sin 2 2(9 rct and 5 MNS . 
They are obtained as follows. We first calculate 
the expected numbers of /j.-like and e-like events 
at HK for the following input parameter values: 





= 3 x 10~ 3 cV 2 , 


(39a) 


6m lo L 


= 7 x 10~ 5 eV 2 , 


(39b) 


sin 2 2(9 ATM 


= 1, 


(39c) 


sin 2 26, sol 


= 0.85, 


(39d) 


sin 2 26> RCT 


= 0.01 or 0.06, 


(39e) 


^MNS 


= 0° , 90° , 180° or 270° , 


(39f) 



by assuming the constant matter density of p = 
3g/(cm) 3 along the baseline. The resulting num- 
ber of events are then analyzed by varying all the 
6 parameters under the known constraints for the 
four parameters 8m 2 , sin 2 29 , , Sm 2 , , and 

1 ATM ' ATM 1 SOL ' 

sin2 20 S oL = 4 ta n 2 9 SOL /(l + ^n 2 # SOL ) 2 , and by 
allowing for 3% normalization uncertainties in the 
earth matter density and the overall beam flux. 
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Figure 5. Allowed regions in the plane of 
sin 2 2t9 RCT and S MNS after 6 years of T-to-HK. 
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In addition, some of the known systematic errors 
from the SK and K2K experiences were taken into 
account. 

We find from Fig. 5 that the T-to-HK experi- 
ment can establish CP violation at the 3-cr level 
if <5 MNS - 90° or 270°. Also, the use of both 2° 
and 3° OAB can contribute to resolving the de- 
generacy between 5 , ~ 0° or 180° if sin 2 29 „ 

° J MNS RCT 

is not too small. When comparing the contours of 
Fig. 5 with the corresponding results of rcf. [20] 
which were obtained for the narrow-band beams 
(NBB), however, we find that the areas of the 3- 
a allowed regions are significantly wider for the 
OAB than the NBB and that the capability of dis- 
tinguishing between the 0° and 180° degeneracy 
is lower for the OAB. These results are obtained 
even though the typical flux of OAB is higher 
than that of NBB, mainly because of much higher 
background of the OAB, especially in the high en- 
ergy tail. The background is found most serious 



for the ' 



measurements where the 



v e com- 



ponent of the 77^ enriched beam contributes to 
the e-like signal since HK cannot distinguish the 
charge. We find that the background level is sig- 
nificant in the high energy tail of the OAB, and 
hence rejection of events with high neutrino ener- 
gies can improve the measurement. Our studies 
on the T-to-HK experiments with OAB will be 
reported in [21] soon. 

4. The parameter degeneracy problem 

Those studies, however, assumed that the 
neutrino mass hierarchy is known to be nor- 
mal (£mf 3 = Sm^ TM ) and also we assumed 
sin 2 8 , = 0.5 as an input. 

ATM L 

In particular, our results for <5 MNS = 90° and 
270° in Fig. 5 suggest that the T-to-HK experi- 
ment can establish CP violation at the 3-cr level 
or higher, if sin 2 2# RCT is not too small. We find, 
however, that this observation depends critically 
on the assumption that we know the neutrino 
mass hierarchy to be normal, i.e., the hierarchy 
pattern I in Fig. 1. 

As a demonstration, we show in Fig. 6(a) the 
allowed region of the T-to-HK simulation by an- 
alyzing exactly the same input 'data', calculated 
for the normal hierarchy (5m 2 3 > 0), by assum- 



(a) 



(8m- M) assumed) 

i , 



(&n? 3 <0 assumed) 



90 180 270 

5mns 



(sin 2 0A™,$in 2 20RCT) C 3 ) 
true mirror 



(045,0.06) 




0.40 0.50 0.60 0.70 

sin 2 ftvrM 



Figure 6. Doubling of the regions allowed by fu- 
ture LBL experiment T-to-HK, between J-PARC 
at Tokai and HyperKamiokande, in the plane 
of S MNS and sin 2 26» RCT 



sin 2 2# RCT (b). 



(a), and sin 9 ATM and 



ing that the neutrino mass hierarchy is inverted 
(cw 2 3 < 0). We find the region encircled by thin 
lines from this analysis. The minimum of the \ 2 
value in those regions are only slightly larger than 
unity, and hence the solutions cannot be excluded 
by the measurements. Most strikingly, the al- 
lowed regions obtained by assuming the inverted 
hierarchy give <5 MNS ~ 0° or 180°, i.e., the very 
conclusion of the evidence of CP violation in the 
lepton sector depends on our capability of dis- 
tinguishing between the two remaining hierarchy 
cases, I (normal) or III (inverted). 
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There is one more hidden degeneracy problem 
among the 6 parameters of the 3 neutrino model. 
It is the magnitude of sin 2 6> ATM = |C/ M 3| 2 which 
has a two-fold ambiguity 

1± J I - sin 2 26* 
-n 2 ATM = M = —^— 2 — (40) 

in the last line of eq. (22). The two- fold ambiguity 
was hidden so far because we chose the extreme 
value of sin 2 20, = 1 as our input. It should 

ATM 1 

be remarked that the presently allowed range of 
sin 2 20 > 0.90 in eq. (1) corresponds to 

0.34 < sin 2 6> < 0.65 . (41) 

ATM v ' 

In Fig. 6(b), we showed the allowed region in the 
same T-to-HK simulation, when the input data 
are calculated for 

sin 2 6> = 0.35 or 0.45 (42) 

ATM v ' 

instead of sin 2 ATM = 0.5 in cq. (39). In addi- 
tion to the correct solutions which arc marked 
by the solid blob at the input points, we find a 
mirror solution at the sin 2 9 , value correspond- 

ATM 1 

ing to the opposite sign in eq. (40) for the same 
sin 2 26* value. We notice also that the mirror 

ATM 

solutions appear at different values of sin 2 20 RCT , 
and hence the measurement of sin 2 29 „ at T2K 

ROT 

and in the future T-to-HK experiments also suffer 
from this degeneracy problem. The reason for the 
correlation between sin 2 # ATM and sin 2 20 RCT is a 
manifestation of the fact that the magnitude of 
the leading term that dictates the — ► v e oscilla- 
tion probability (14c) is their product in eq. (17). 
This correlation suggests that if we can measure 
sin 2 29 RCT that governs the leading term of the 
v e — ► v e oscillation in cq. (14b) accurately inde- 
pendent of the magnitude of sin 2 9 ATM , we can 
resolve the ambiguity. 

These studies demonstrate clearly that we need 
more works to do in order to determine the neu- 
trino mass hierarchy pattern between the normal 
(I) and inverted (III) ones, and to determine the 

third mixing angle sin 2 29 and sin 2 9 inde- 
ed G RCT ATM 

pendently. New generations of the reactor neu- 
trino experiments, such as Double-CHOOZ [22] 
and KASKA [23] have been proposed to tackle the 



latter problem. In the rest of this talk, I would 
like to introduce our studies about possible future 
long baseline neutrino oscillation experiments in 
Asia, which make use of neutrino beams from the 
J-PARC facility. 

5. Fate of the T2K and T-to-HK beam 

We first study the fate of the off axis beam of 
the T2K and the possible T-to-HK project. The 
OAB at J-PARC is shoot underground, such that 
exactly the same OAB will be available at the 
present SK site and the planned HK site. The 
center of the beam at around Kamioka site passes 
through more than 10 km deep in the earth crust. 
The beam center will re-appear on the earth sur- 
face in the middle of the Japan sea, and the OAB 
of the off-axis angle between 0.5° and 3° will be 
available free in Korea. The surface view of the 
OAB (2.5°) is shown in Fig. 8 and the vertical 
cross section view of the beam is shown in Fig. 9. 

In the surface view of Fig. 8, the red curves 
show the contours of the same baseline length. 
The solid and dashed black contours give the off- 
axis angle. The center of the beam is at 0°, which 
appears in the Japan few hundred km east 

from the Korean coast. Inside of Korea main- 
land, the off-axis beam of 1° to 3° is available at 
L =1,000 km to 1,200 km. The off axis beam at 
0.5° is available at the Korean east coast, for the 
OAB(3°) at Kamioka. 

Therefore, depending on the location of a pos- 
sible neutrino detector in Korea, the flux of off- 
axis angle between 0.5° and 3° shown in Fig. 4 
will be available during the period of the T2K 
experiment, and possibly during the T-to-HK ex- 
periment. Because the flux at L <~ 1,000 km to 
1,200 km is about a factor of 10 to 16 smaller 
than that at Kamioka, respectively, we may need 
a detector of 100 kton volume in order to obtain 
significant new contributions. 

Because the OAB at large angles has a sharp 
peak at low energies, see Fig. 4, one may study 
the oscillation probabilities at large L/E. There 
are a few advantages of neutrino oscillation exper- 
iments at large L/E while keeping the neutrino 
energy at sub-GeV. 

• A12 grows with L/E, and hence the sensi- 
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tivity to <5 MNS in the — » ^ e experiment 
will grow; see eq. (14c). 

• The sign of the term linear in both A12 and 
A13 determines the neutrino mass hierar- 
chy, because the term A12 • A13 is positive 
for the normal hierarchy (I), while it is neg- 
ative for the inverted hierarchy (III). 

Relative smallness of the earth matter effects at 
lower energies will be helpful in the study since 
the above features of the — > v e oscillation prob- 
ability can be read off from the expansion in terms 
of A12 in eq. (14c) which is a good approximation 
only up to a few GeV; see Fig. 2. 

In Fig. 7, we show the oscillation probabilities 
P{v» -» fe) and P(V^ -» 77 e ) for E v = 0.7 GeV 
plotted against the baseline length L. The up- 
per two figures are for the normal hierarchy (I) 
and the lower two are for the inverted hierarchy 
(III). We can see from the figure that in some 
cases, small differences at L = 300 km for the 
T2K and T-to-HK experiments can be enhanced 
significantly in Korea between L =1,000 km and 
1,200 km. The difference between the normal 
and inverted hierarchy can be very significant and 
there is a hope that a combination of the T-to- 
HK measurement and the T-to-Korea experiment 
can resolve the ambiguity in the <5 MNS and the neu- 
trino mass hierarchy cases at the same time. The 
study is underway [24] . 

At the moment, we do not know if there is a 
strong interest in constructing a huge neutrino 
detector in Korea. Because we feel that this in- 
vestigation is rather promising, we study the fate 
of the OAB from J-PARC in the T2K experi- 
ment carefully. Although only the fate of the 
OAB(2.5°) is shown in Fig. 8 and Fig. 9, we made 
similar figures for OAB (2°) and OAB (3°). In case 
of OAB (2°), the smallest off- axis angle at the east 
coast of Korea is about 1.5°, while for OAB(3°), 
the smallest angle is about 0.5°. 

In the cross section view of Fig. 9, we show the 
baseline, along the direction of J-PARC at Tokai 
village to the center of the SK and the proposed 
HK site in Kamioka. The OAB(2.5°) beam hits 
SK and HK at the same off-axis angle of 2.5°, 
while the beam center (0°) reaches the earth sur- 
face in the sea. It is rather important to rec- 
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Figure 7. The oscillation probabilities P{v^ — > 
v e ) and P(V^ — » V e ) for E v = 0.7 GeV plotted 
against the baseline length L. The upper two 
figures arc for the normal hierarchy (I) and the 
lower two are for the inverted hierarchy (III). 
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Figure 8. The fate of the OAB 2.5 degree beam from J-PARC. Surface view. 




Figure 9. The fate of the OAB 2.5 degree beam 
from J-PARC. Vertical view. 



ognize that the Japan sea is rather shallow, and 
the beams which reach Korea do not travel un- 
der water. The numerics in the earth crust show 
the effective mass density p in units of g/(cm) 3 . 
We can view that under the sea the earth crust is 
thin, and the mantle appears already at 20 km be- 
neath the earth surface. The density variations, 
however, are not expected to be significant. We 
hope that our studies [24] will be found useful by 
our colleagues in Korea. 



6. High energy super beam to China 

Although we do not yet know if there are strong 
enough interests in constructing a huge neutrino 
detector in Korea, strong interests have been ex- 
pressed by our Chinese colleagues about the pos- 
sibility of sending super neutrino beams from J- 
PARC at Tokai to somewhere in mainland China. 
A possible 100 kton level water Cerenkov detector 
BAND (Beijing Astrophysics and Neutrino De- 
tector) [25] has been proposed, and if it will be 
placed in Beijing, the baseline length from Tokai 
will be about L=2,100 km. The unique capability 
of the BAND detector is that it is a segmented 
detector, where each segment contains 10 ton of 
water surrounded by the photo-multiplier tube to 
measure the Cerenkov light, and at the same time 
has a good capability of calorimetric energy mea- 
surement [25]. 

Encouraged by the interests expressed by our 
Chinese colleagues, we made a serious study of 
physics capability of a very long baseline neu- 
trino oscillation experiments between J-PARC at 
Tokai and BAND in Beijing [11]. The global view 
of the eastern Asia seen from J-PARC at Tokai 
is shown in Fig. 10. It is clearly seen from the 
figure that the LBL experiments with L greater 
than 2,000 km are feasible if a huge neutrino de- 
tector can be built in China. We can also tell 
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Figure 10. East Asia viewed from J-PARC 
in Tokai. The baseline length L = 295 km 
(Kamioka), 1,200 km (Seoul) and 2,100 km (Bei- 
jing). 



from the figure that the direction from Tokai to 
Kamioka is almost exactly from east to west, and 
the T2K beam goes through southern part of Ko- 
rea. Since Beijing is placed north east from the 
T2K beam-line, a new neutrino beam line toward 
Beijing should be constructed at J-PARC. 

In addition, in order to send sufficient neutrino 
flux to the distance greater than 2,000 km, and 
in order to study physics at the same L/E region 
where the first oscillation mode takes its maxi- 
mum value, 

|A 13 | ~ 7T when L(km)/E(GeV) ~ 500, (43) 

for 5m 2 , = 2.5 x 10~ 3 eV 2 , we need a neutrino 

ATM ' 

beam whose flux peaks at around 3 to 6 GeV. 
Such neutrino beams can indeed be designed for 
the J-PARC 50 GeV proton accelerator, and has 
been named HENBB (High Energy Narrow Band 
Beam) in ref. [11]. We show in Fig. 11 the ex- 
pected flux times the neutrino energy E v as func- 
tions of E v . We multiply the flux by E v in or- 
der to show the peak in the expected number of 
events in the absence of the neutrino oscillation. 
The peak position of the flux xE tJ is called E p 



Figure 11. High energy narrow band beams 
(HENBB) at J-PARC for the possible T2B 
(Tokai-to-Beijing) project. 



and the simulation results are shown for E p = 
3.3, 4.3, 5.7 and 6.4 GeV. The optics adopted for 
generating the HENBB makes the peak rather 
sharp, especially in the high energy side. This 
sharp cut-off of the high energy neutrinos makes 
the background level low when the HENBB is 
used for the LBL experiments. In fact we find 
that the largest background from the high energy 
neutrinos beyond E p comes mainly from the con- 
tribution from the K — > [iv decays, which may be 
measured precisely at a near detector. We made 
a parameterization of the HENBB flux for a given 
E p in order to save time for the beam simulation, 
and the parameterization is given in ref. [11]. 

We show in Fig. 12, the expected total numbers 
of the e-like events and the //-like events. Since 
we do not require the charge identification capa- 
bility for the BAND detector* the //-like events 
are the sum of and charged current (CC) 
events and the e-like events are the sum of v e 
and T> e CC events and those expected from the 
NC events where produced ir°'s are mistaken as 
e^. At higher energies, t production from the 

t Possibility of putting the segmented water Cerenkov de- 
tector of BAND under the magnetic field has been studied 
by the IHEP group [25]. 
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(a) NBB (E P eak=4GeV) at L=2100 km 
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high energies, because such significant suppres- 
sion has not been identified in our previous stud- 
ies for the T-to-HK project, where the obser- 



10 20 30 40 50 60 70 80 

N(e)/500kton-year 



Figure 12. The expected number of mu-like and 
e-like events for the future T2B experiment with 
the beam HENBB(4GcV). 



dominant —* v T oscillation mode becomes sig- 
nificant, and their leptonic decay modes have also 
been counted as background. The contours shown 
in Fig. 12 are the prediction of the three neutrino 
model for 5m 2 , = 3.5 x 10 -3 eV 2 , which was the 

ATM ' 

preferred value a few years ago. The results are 
for 500 kton-year, or for 5 years with the 0.75 MW 
J-PARC and a 100 kton BAND. They are shown 
for all the hierarchy cases I, II, III, IV of Fig. 1, 
where the results for the cases II and IV can be 
interpreted as the expectation for anti-neutrino 
oscillation experiments, according to the theorem 
eq. (11), after the ratio o CG {pi>)lo cc (vi) ~ 
is multiplied,. 

It is most striking to learn from the figure 
that the — > v e oscillation probability is al- 
most a factor of 7 smaller in the inverted hi- 
erarchy (III) than the normal hierarchy (I). On 
the contrary, the — > T7 e oscillation is sup- 
pressed for the normal hierarchy (IV). This is 
clearly due to the enhanced matter effects at 



vation of the v u 



v e events arc essential to 



determine the CP-violating parameter, J MNS or 
sin 5 M • The two experiments measure approx- 
imately the same phase of the oscillation proba- 
bility, because L/E = 295 km/0.7 GeV is not far 
from 2,100 km/4 GeV. 

This can be understood qualitatively as follows. 
The oscillation probability P VjM ^ Ve inside of the 
matter can be expressed as eq. (14c), where the 
mass squared differences and the mixing matrix 
are replaced by the corresponding ones Sfh 2 j and 

U a i according to eq. (23). From eq. (14c), we 
can tell that the difference in P v ^ Ve between 
the mass hierarchy I and III should come from 
the difference in the second term proportional to 
A12, which changes the sign between I and III. We 
find that the mixing matrix terms do not change 
significantly in the relevant energy region, and 
the second term contributes constructively for the 
normal hierarchy (I), whereas destructively for 
the inverted case (III), as long as | Ai 3 | < ir. From 
Fig. 2, we find that the magnitude of the term 
8fh\ 2 grows rapidly in a few to several GeV re- 
gion, approaching that of the leading terms of the 
order Sm 2 . Although the similar trend is ob- 

ATM G 

served for the oscillation at E v ~ 0.7 GeV for the 
T2K and T-to-HK projects, the effects are small 
enough to be compensated by shifts in the other 
model parameters, notably sin 2 2# RCT and 5 MNS - 
This drastic suppression of — ► v e oscillation 
probability for the inverted hierarchy at high en- 
ergies gives us a powerful tool to determine the 
neutrino mass hierarchy at very long baseline neu- 
trino oscillation experiments which are capable of 
detecting electrons. The difference in the predic- 
tions of the normal and the inverted hierarchy 
cases at sin 2 2$ ~ 0.1 is so huge that no ad- 
justment of the other model parameters can com- 
pensate for it. The Tokai-to-Beijing project with 
HENBB seems to be an ideal example of such ex- 
periments. In Fig. 13, we show the results of sim- 
ple statistical analysis [11], where we calculated 
the number of expected events at BAND after 5 
years of HENBB at E p = 4 GeV and another 5 
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Hierarchy I (true) vs III (fit) 




0.02 0.04 0.06 0.08 0.1 0.12 



Figure 13. The expected minimum chi-squared of 
the combined T2K and T2B experiments when 
the normal hierarchy is realized in the nature 
while the inverted hierarchy is wrongly assumed 
in the analysis. 



years at 6 GeV for the normal hierarchy and then 
by analyzing the 'data' by assuming the inverted 
hierarchy. The results of our analysis depend se- 
riously on the assumed value of sin 2 20 R , and 
hence the Xmin vames °f our statistical analysis 
is shown as functions of the sin 2 29 „ value in 

RCT 

the fit. We examined the cases with the input 
('true') parameters at sin 2 20£™ = 0.1, 0.04 and 
0.02, for the 4 phase angles, 4™ = 0°,90°, 180° 
and 270°. It can be seen from the figure that 
if sin 2 20^™ > 0.04 and the hierarchy is nor- 
mal (I), then the inverted hierarchy case (III) 
can be excluded at 3-cr level for any values of 
<S TC . On the other hand, the generated 'data' 

MNS ' ° 

at sin 2 20*,™ = 0.02 for ^™ = 0° and 90° can be 
fitted well within the allowed range of the other 
model parameters. 

We also studied the possibility of constrain- 
ing the S M parameter in the possible Tokai-to- 
Bcijing project by using the HENBB, but the re- 
sults we find is not as promising as those from the 
analysis [20] of the possible Tokai-to-HK project. 



The reason is rather simple. If the neutrino mass 
hierarchy is normal, the oscillation probability for 
the 77^ — > 77 e is strongly suppressed, and hence we 
cannot expect much number of events by sending 
anti-neutrino beams. If the hierarchy is inverted, 
on the other hand, the — > v e oscillation prob- 
ability is strongly suppressed, and we will need 
significantly higher power neutrino beams. If the 
inverted hierarchy is realized in the nature, the 
Tokai-to-Bcijing project with a 4 MW proton ac- 
celerator at J-PARC may be necessary. In the 
Tokai-to-HK project, the use of the low energy 
neutrino beam helps reducing of the matter ef- 
fects and a one Mton detector more than compen- 
sates for the small cross section at low energies. 

In rcf. [11], we show the results of our analysis 
at L = 1, 200km with the HENBB at E p = 3 GeV 
and 5 GeV, and find nearly as promising results 
as the above Tokai-to-Bcijing analysis. If there 
appears a strong interest in constructing a huge 
neutrino detector in Korea, this might as well be 
an appealing option. There, the remnant of the 
OAB will be studied during the period of the T2K 
and the planned T-to-HK project, and at a later 
stage, a new beam-line for the HENBB may be 
constructed at the J-PARC site to determine the 
neutrino mass hierarchy. 

7. Neutrino Factory at J-PARC 

In the final section of my report, I would like 
to introduce our studies on the possibility of the 
future very long baseline (VLBL) neutrino oscil- 
lation experiments with a neutrino factory at J- 
PARC and a BAND like 100 kton level detector 
at a few thousand km away. The possibility of 
a neutrino factory at J-PARC is being studied 
in Ref. [26]. For definiteness, we assume a 100 
kton-lcvcl segmented water- Cerenkov calorimeter 
detector BAND at L = 2, 100km away from J- 
PARC, and a neutrino factory at J-PARC which 
is capable of delivering 10 21 n + or decays at 
10 GeV in one year. We will show that the goals 
of determining all the neutrino model parameters, 
including the neutrino mass hierarchy, the sign of 
sin 2 ATM — 1/2, and the degeneracy in the CP 
phase 5 MNS , can be reached if such an experiment 
can be realized. 
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The physics prospects of VLBL oscillation ex- 
periments with a neutrino factory has been stud- 
ied in the past by assuming that the detector can 
identify charges [27], and hence the possibility of 
charge identification at BAND has been inves- 
tigated [25]. We have tried to show in the re- 
port [28] that even if the detector is charge blind, 
the VLBL experiment will achieve all the physics 
goals if it is capable of distinguishing the /i-like 
and e-like events with high confidence level, and 
if the detector is capable of determining the event 
energy (the incoming neutrino energy minus the 
outgoing neutrino energies) at the 10 % level. In 
the following, we assume that the proposed de- 
tector BAND will have such capabilities. 

In the neutrino factory, neutrinos are produced 
from the decay of high energy muons, /i+ — > 
v^v^ or /i~ — > V[J; e e~ . Not only the shape 
and the ratios of the neutrino fluxes but also their 
overall normalization will be known accurately in 
the neutrino factory experiments. Assuming the 
relativistic muons, the 17^ and v e (y^ and V e ) 
fluxes from /x + (/i~) beam are expressed as 



ttL 2 
7 2 n 



2y 2 [(3-2 2/ ) T ^(l-2 2 /)],(44a) 



*„.c*.) = -^£1 V[(i -v)t W - v)\ > ( 44b ) 

where 7 = E^/m^, y = E„/E^ 7 and P M is the 
average muon polarization, and is the number 
of the decaying muons for which we assume 10 21 
per year in the following. 

The e-like signal from /i + beam, A r e (/i + ), is 
given by the sum of e + from — ► V e appear- 
ance and e~ from the v e — ► v e survival mode, 
whereas the /U-like signal, A^(/i + ), is the sum of 
/i + from — > and jjT from v^—^v^ 

N e {ji + ) : (u e -» v e ) + (77 M -> 77 e ) , (45a) 
JV M (/i+) : (*V -> 77 M ) + (i/ e -» . (45b) 

The signals from the /U~ beam, N e (fj,~) and 
Nn(fj,~), are obtained in the same way; 

N e {n~) : (V e -> 77 e ) + (^ -> i/ e ) , (46a) 

iV M (/x-) : -» + (V e > F^) . (46b) 

These signals in the i-th energy bin, 
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(47) 



are then calculated by ignoring the energy res- 
olution errors, but correcting for the missing 
neutrino energies in the events where and 
comes from t ± production via the charged 
currents. In the following analysis we dropped 
from our analysis low energy events with E b s < 
2 GeV. 

The results are shown in Fig. 14 for 10 21 de- 
caying /i + (top two figures) and n~ (bottom 
two figures) each at 10 GeV. The results for 
the normal-hierarchy are given in thick lines and 
those of the inverted hierarchy arc in thin lines; 
the solid, dashed, and dot-dashed lines are for 
sin 2 ATM =0.5, 0.35, and 0.65, respectively. The 
other model parameters are chosen as in cq. (39) 
atsin 2 2fl RCT =0.06and S MNS = 0°. 

The contribution from the survival mode dom- 
inates the /n-like signal but they vanish at around 
E„ ~ 5 GeV for both A^(m + ) and A^(m~), where 
most of [i neutrinos oscillate into r neutrinos. 
Therefore the /i-like events in the surrounding en- 
ergy bins are sensitive to the v e — ► or 77 e — > 
transition probability. 

For instance, in the histogram for N fl (fj, + ), we 
can see significant dependence on the sign of 
sin 2 # ATM — 1/2, especially around E v ~ 5 GeV. 
The sensitivity is due to the v e — > appear- 
ance mode, whose probability is proportional to 
sin 2 ATM sin 2 # RCT , just like the — ► v e oscilla- 
tion probability in eq. (14c). In order to resolve 
the degeneracy in the sign of sin 2 — 1/2, how- 
ever, we need to have an independent measure- 
ment of sin 2 26> RCT . 

We observe from Fig. 14 that N e (/j, + ) is very 
insensitive to sin 2 ATM — 1/2 throughout the en- 
ergy bins. This is especially so for the normal hi- 
erarchy, where the — > V e oscillation is strongly 
suppressed due to the matter effects, as explained 
in the previous section, see Fig. 12. Although the 
suppression is absent for the inverted hierarchy, 
the contributions from the F M — > 77 e transition is 
still suppressed significantly because of the small 
anti- neutrino cross sections on the matter target. 
It is therefore envisaged that the measurement of 
N e (fj, + ) will give us the desired independent mea- 
surement of sin 2 20 r:rr , because it is dominated by 

RCT J 

the v e — > v e survival mode, especially when the 
mass hierarchy is normal. Therefore, the neutrino 
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factory VLBL experiments with a detector capa- 
ble for distinguishing events from ^ can give 
us a precise measurement of sin 2 2# RCT just like 
the next generation reactor experiments [22,23]. 
In addition, the presence of the enhanced large 
matter effects for the <-> v e and v y, <-> F e make 
the inter-dependences of the four observables in 
Fig. 14 very different between the normal and 
the inverted hierarchies. We also find that the 
results depend rather significantly on the input 



value if sin 

MNS 



20 T is not too small. Because 



2 3 4 5 



7 8 9 10 



the energy dependence of the signal can be stud- 
ied in the wide range of L/E at VLBL exper- 
iments with a neutrino factory, the degeneracy 
between 5 =0° and 180° can also be resolved 

MNS 

[28]. 

We performed a \ 2 analysis similar to the ones 
we performed for the Tokai-to-HK [20] and the 
Tokai-to-Beijing [11] projects, in order to study 
these questions quantitatively. The additional as- 
sumption that we make for the detector is that it 
is capable of measuring the event energy calori- 
mctrically, with the accuracy which makes the er- 
ror of our analysis based on the 1 GeV bin his- 
togram small. The energy threshold of E v > 
2 GeV has been introduced to avoid numerical 
sensitivity to the binning due to rapid oscillation. 
Our analysis is hence performed for 8 bins be- 
tween E v = 2 GeV and 10 GeV. 

Because the expected number of events is huge 
for the above experimental set up, a few tenths of 
sounds of events for each bin, as can be seen from 
Fig. 14, we introduced the following systematic 
uncertainties in the analysis: 

• 2 % error each in the detection efficiency of 
e-like and /Lt-like events 

• 2 % error each in the CC cross section of 
neutrinos and anti-neutrinos 



Figure 14. Number of e- and /^-like signals for 
10 21 decaying [i + (top two) and \i~ (bottom two) 
each at 10 GeV. Thick lines are for the nor- 
mal hierarchy, and thin lines are for inverted; 
sin 2 9 ATM =0.5 (solid lines), 0.35 (dashed lines), 
and 0.65 (dot-dashed lines). 



Here we assumed that the errors in the detection 
efficiencies are independent for e-like and /z-like 
events, and also the errors in the v and V cross 
sections are independent. On the other hand we 
assume that the errors for v e and cross sec- 
tions are common (100 % correlated), as well as 
for their anti-neutrino counterparts. Since the ex- 
periments will be performed more than a decade 
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Figure 15. Xmin as a function of the fitting param- 
eter sin 2 20 RCT when the input events are gener- 
ated for the normal hierarchy (I) and the analysis 
is done by assuming the inverted hierarchy (III). 
The input data are calculated for sin 2 A m ° =0.5 
at sin 2 20*™ = 0.1, 0.008 and 0. The four curves 
are for 5 MNS = 0° (solid), 90° (dotted), 180° 
(dash) and 270° (dot-dash). 
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Figure 16. Xmin as a function of the fitting 
sin 2 ATM • The input data of the events are cal- 
culated for sin 2 0*^™ =0.35, 0.4, and 0.45 with 
sin 2 26^ c u ° = 0.06 (solid lines), 0.08 (dotted lines) 
and 0.1 (dashed lines) and <^™ = 0°. The other 
input values are the same as in eq. (39). 



in the future, it is possible that the systematic er- 
rors can be reduced even further. In ref. [28], we 
also show the results when the above systematic 
errors are set to zero. 

We show in Fig. 15 the X m in vame as a function 
of sin 2 20 RCT , when the 'data' are generated by 
assuming the normal hierarchy, while the anal- 
ysis is done by assuming the inverted hierarchy. 
We show our results for sin 2 20^™ = 0.1, 0.008, 
and 0, and for the four CP violating phase an- 
gles; 5 MNS = 0°, 90°, 180° and 270°. The results 
show that we can distinguish between the normal 
and the inverted mass hierarchy at 3-ct level, if 
sin 2 20*™! > 0.008. At small sin 2 26* nrr , we find 
that the contributions of the transition modes be- 
come negligible, but the e-like signals from both 
/i + and fj,~ beams are useful in determine the hi- 
erarchy. This is because the disappearance prob- 
ability is suppressed for the V e — ► V e transition in 
the normal hierarchy, whereas that of v e — > v e 
is suppressed in the inverted hierarchy. How- 
ever, the differences of the e-like signals between 



the normal and the inverted hierarchy reduces to 
<~ 1% when sin 2 20 RCT = 0.004, and the sensi- 
tivity is limited by our knowledge on the ratio 



'CC 



>e)/o- C cOe)- 



We show in Fig. 16 the X m in as a function of 



sin 2 61 



The input 'data' are calculated for 



sin 2 6f T u ° = 0.35, 0.4, and 0.45 for three values 
of sin 2 20 R r ™ = 0.1 (dashed lines), 0.08 (dotted 
lines) and 0.06 (solid lines), and for = 0° 

with the normal hierarchy. The values of the 
other parameters are taken as in cq. (39). The 
X^^function is found by varying the fitting pa- 
rameters within the normal hierarchy. We see 
that each x^has two dips at the sin 2 atm values 
which give the same sin 2 20 ATM . The results show 
that we can resolve the degeneracy at the 3-cr level 
when sin 2 0' rue = 0.35 (sin 2 20 A = 0.91) or 0.4 

ATM V ATM ' 

(sin 2 20 ATM = 0.96) for all the three sin 2 20*™° 
values considered (sin 2 20 R "° > 0.06), but that 

truc = 0.45 (sin 2 
resolved from sin 2 9 , 



sin 2 



sin 2 20„ 



= 0.99) can be 
0.55 at 3-cr only if 



>0.1. 



In Fig. 17, we show the allowed regions in the 
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Figurc 17. Allowed regions in the plane of 
sin 2 20 RCT and 5 MNS for sin 2 26> R r c u ° = 0.02, 0.06 
and 0.1, for 5£™ = 0°, 90°, 180° and 270°. 



_,2 o^truc 

RCT 

0.02 , 0.06 , and 0.1 and 5£™ = 0° (top), 90° (sec- 
ond), 180° (third), and 270° (bottom) in the nor- 
mal hierarchy. In each figure, the input parame- 
ter points are shown by solid-circles. The normal 
hierarchy is assumed in the fitting. The regions 
where Xmin< 1; 4, and 9 are depicted by solid, 
dashed, and dotted boundaries, respectively. 
The figures for S true = 90° and S true = 270° 

° MNS MNS 

show that the CP phase 5 MNS can be constrained 
locally around the 'true' points for all the input 
sin 2 2# RCT values. We can hence establish CP vio- 
lation in the leptonic sector, and discriminate the 
maximal CP violation, 5 l ™° = 90° or 270°, from 
the CP conserving cases of <5 MNS = 0° or 180° at 
the 3-er level if sin 2 20 R r c u ° > 0.02. The discrimina- 
tion between <5 MNS = 0° and 180° is also possible 
at the same level, but for the small 4 < x 2 nin < 9 
island at sin 2 2(9*^ = 0.02. 

The above results for CP violation are compa- 
rable to those of Tokai-to-HK studies shown in 
Fig. 5. The sensitivity in this analysis, with the 
neutrino factory and a charge-blind detector, is 
based essentially on the T violating difference be- 
tween P v ^ V)i and P„ _n, e . Because the matter 
effects enhance both of these transitions whereas 
suppress strongly the transitions — > v e and 
v e — > Vp, N fl (fi + ) is sensitive to P Ve ^> Vll and 
N e (fi~) is sensitive to P v ^ Vc - Therefore, the 
above results are obtained only when the neutrino 
mass hierarchy is normal. A significantly less ac- 
curate results are expected if the mass hierarchy 
turns out to be inverted, because the enhanced 
transition modes appear for the anti-neutrinos 
whose cross sections on the matter is a factor of 
3 smaller than the neutrino cross sections. 

Summing up, the following results are expected 
for a VLBL experiment with a neutrino factory 
which delivers 10 21 decaying [i + and /i~ at 10 
GeV and a 100 kton detector which is placed 
2,100 km away and is capable of measuring the 
event energy and distinguishing from fi^, 
but not their charges: The neutrino mass hier- 
archy can be determined if sin 2 



the degeneracy in sin 



26 RCT > 0.008, 
can be resolved for 



sin 2 2# ATM = 0.96 if sin 2 20 RCT >0.06, and the 

can be uniquely con- 



ATM 

CP-violating phase 5 
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strained for sin 2 29 r , > 0.02 if its true value is 

RCT <~ 

around 90° or 270°, while it can be constrained 
for sin 2 2# RCT > 0.04 if its true value is around 0° 
or 180°, all at the 3-cr level. 
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